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We investigate correlations of the Polyakov loop fluctuations with eigenmodes of the lattice Dirac 
operator. Their analytic relations are derived on the temporally odd-number size lattice with the 
normal non-twisted periodic boundary condition for the link-variables. We find that the low-lying 
Dirac modes yield negligible contributions to the Polyakov loop fluctuations. This property is 
confirmed to be valid in confined and deconfined phase by numerical simulations in SU (3) quenched 
QCD. These results indicate that there is no direct, one-to-one correspondence between confinement 
and chiral symmetry breaking in QCD in the context of different properties of the Polyakov loop 
fluctuation ratios. 


I. INTRODUCTION 

Color confinement and chiral symmetry breaking are 
the striking non-perturbative phenomena in low-energy 
QCD, which are of particular importance in particle and 
nuclear physics [ll-Q. 

Several scenarios of the confinement mechanism have 
been proposed , where the ghost and gluon propaga¬ 
tors in deep infrared need to be quantified, thus this re¬ 
quires a non-perturbative quantization of QCD. Whereas 
this issue has been investigated extensively, it remains 
challenging to comprehend the non-perturbative aspects 
from the first-principle calculations. 

In a pure SU(3) gauge theory, the Polyakov loop is 
an exact order parameter of the Z 3 center symmetry 
and for deconfinement, which dictates a first order phase 
transition SiOl- In the presence of light dynamical 
quarks, the Polyakov loop loses its interpretation as an 
order parameter and is smoothly changing with temper¬ 
ature. However, contrary to the broad Polyakov loop, 
a particular ratio of the Polyakov loop susceptibilities 
retain a clear remnant of the underlying Z 3 center sym¬ 
metry fairly well even in full QCD with the physical pion 
mass [13, [3 • Thus, the ratio of the Polyakov loop fluc¬ 
tuations can serve as observables to identify the onset of 
deconfinement in QCD. 

In the presence of light dynamical quarks, the tran¬ 
sition from hadronic phase to quark-gluon plasma be¬ 
comes crossover and accompanies partial restoration of 
chiral symmetry at finite temperature [iH3. Spon¬ 
taneous chiral symmetry breaking is characterized by a 
non-vanishing condensation of quark-bilinear operators. 
The low-lying Dirac modes, which are the eigenmodes 
of the Dirac operator with small eigenvalues, are known 
to be responsible for saturating the chiral condensate 
of li ght quarks through the Banks-Casher rela¬ 

tion [T3|. 

In fact, at vanishing and small baryon chemical poten¬ 
tial, the lattice QCD results suggest that there is an inter¬ 
play between quark deconfinement and chiral crossovers 
as they take place in the same narrow temperature range 


[l3| . Also, in the maximally Abelian gauge, confinement 
and chiral symmetry breaking, as well as the instantons, 
simultaneously disappear when the QCD monopoles are 
removed |19l-l^. On the other hand, there exist sev¬ 
eral observations that unbroken chiral symmetry does 
not dictate deconfinement: Given a tower of hadron spec¬ 
tra with eliminating the low-lying Dirac modes [^ . the 
hadrons keep their identity even in chirally restored phase 
where parity doublers are all degenerate. In addition, it 
has been shown in the SU(3) lattice simulations that the 
low-lying modes have little contribution to the Polyakov 
loop and to the confining force, indicating that the two 
features are rather independent [23 - 1^ . 

In the context of the interplay between confinement 
and chiral symmetry breaking [I^ [13, [13, [HI, I23l - l33| , it 
is important to make a reliable separation of one from an¬ 
other, whereas those phenomena are supposed to be cor¬ 
related. The apparent coincidence in change of proper¬ 
ties of the Polyakov loop fluctuation ratios and the chiral 
condensate and its susceptibility near the chiral crossover 
might indicate that there is a tied relation between con¬ 
finement and chiral symmetry breaking in QCD. How¬ 
ever, such a possible relation has not been quantified yet. 

Utilizing the Dirac-mode expansion method formu¬ 
lated on the lattice [IJ, the low-lying modes can be sys¬ 
tematically removed in calculating expectation values of 
different operators. 

In this paper, we apply the above expansion method to 
investigate the relation between confinement and chiral 
symmetry breaking in terms of the Polyakov loop fluctu¬ 
ations and their ratios. We put our particular attention 
to the contribution of the low-lying Dirac modes to the 
Polyakov loop fluctuations. 

We derive the analytic relations between the real, 
imaginary and modulus of the Polyakov loop and their 
fluctuations with the Dirac modes on the temporally odd- 
number size lattice with periodic boundary conditions. 
These analytical relations are applicable to both full and 
quenched QCD. We show, through numerical simulations 
on the lattice in quenched QCD, that the low-lying Dirac 
modes yield negligible contribution to the Polyakov loop 
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FIG. 1: The temperature dependence of the Polyakov loop 
susceptibilities ratio Ra = XaIxt from Eq. (jd]) and the chi¬ 
ral condensate of light quarks normalized to its zero 

temperature value. The lattice QCD Monte Carlo results are 
from Refs. [T^ and [^, respectively. The horizontal dashed 
lines are the limiting values of Ra in a pure SU(3) lattice 
gauge theory [T^. The vertical dashed lines indicate the chi¬ 
ral crossover temperature and the range of its errors . 


fluctuations. With these results, also not observed is a 
direct relation between confinement and chiral symmetry 
breaking in QCD through the Polyakov loop fluctuation 
ratios. 

The paper is organized as follows: In the next section 
we derive a set of analytic relations linking the Polyakov 
loop fluctuations to the Dirac eigenmodes. In Sec. Ill, 
we examine the role of the low-lying Dirac modes in the 
Polyakov loop fluctuations and present our numerical re¬ 
sults within quenched lattice QCD. Section IV is devoted 
to summary and conclusions. 


II. POLYAKOV LOOP FLUCTUATIONS 


We utilize the SU(7Vc) lattice QCD formalism and con¬ 
sider a square lattice with spacing a. Each site is indi¬ 
cated by s = (si, S 2 , S 3 , S 4 ) with s^ = ,iV^. A 

gauge held, A^(s) G SU{Nc) with the gauge coupling g, 
is introduced as a link-variable, U^{s) = We 

use spatially symmetric lattice, i.e., Ni = N 2 = N 3 = No- 
and N 4 = AC, with AC > Nr- 

For each gauge configuration, the Polyakov loop L is 
defined as 


1 1\ T ± 

^ = n ^4(s + i4)}, (1) 

s i =0 


where fl is the unit vector in direction of /r in the lattice 
unit, and V is the 4-dimensional lattice volume, V = 


N^Nr. 


Under the Z 3 rotation, the Polyakov loop is trans¬ 
formed into 


L = (2) 

with k = 0, ±1 [ 13 , [ 13 . In confined phase, fc = 0 is 
taken. In deconfined phase, k is chosen such that the 
transformed Polyakov loop L lies in its real sector. 

We introduce the Polyakov loop susceptibilities. 


AfS 

T\A = j^[{m-{\L\n 

(3) 

T^xL = j^[{{LLr)-{LLn 

(4) 

T^Xt = ^[((T^t)^) - {Lt)% 

(5) 


where = Re(L) and Lt = Im(L), and consider their 
ratios, 

Ra = (6) 

XL 

Rt = —. (7) 

XL 


The Polyakov loop susceptibility ratios dH) and 0 were 
shown to be excellent probes of the deconfinement phase 
transition in a pure gauge theory [Hill. They are al¬ 
most temperature independent above and below the tran¬ 
sition and exhibit a discontinuity at the transition tem¬ 
perature. This characteristic behavior is understood in 
terms of the global symmetry of the Yang-Mills La- 
grangian and the general properties of the Polyakov loop 
probability distribution | 12 | . 

In the presence of dynamical quarks, the Polyakov loop 
is no longer an order parameter and stays finite even in 
the low temperature phase. Consequently, the ratios of 
the Polyakov loop susceptibilities are modified due to ex¬ 
plicit breaking of the Z 3 center symmetry. Indeed, both 
Ra and Rt vary continuously with temperature across 
the chiral crossover, however Ra interpolates between the 
two limiting values set by the pure gauge theory. This 
property of Ra is illustrated in Fig. [T] Also seen in this 
figure is that, in spite of smoothening effects observed 
in the presence of quarks, there is an abrupt rate change 
with T in Ra near the chiral crossover T ~ 155 MeV 
The ratio Ra in Fig. [T]has its inflection point at T ~ 150 
MeV, which is fairly in good agreement with the chiral 
crossover range calculated in lattice QCD. Such behavior 
of the Polyakov loop fluctuation ratio can be regarded as 
an effective observable indicating deconfinement proper¬ 
ties in QCD (T^ . 

The apparent modification of Ra and Rt near chi¬ 
ral crossover may suggest that there are certain correla¬ 
tions between confinement and chiral symmetry break¬ 
ing. Such correlations can be best verified when expand¬ 
ing the Polyakov loop and its fluctuations in terms of the 
Dirac eigenmodes. 
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In the following, we formulate the relevant quantities, 
based on this expansion method, and study the influence 
of the low-lying Dirac modes on the properties of the 
Polyakov loop fluctuation ratios. 


III. DIRAC-MODE EXPANSION 


To derive the analytic relation between the Polyakov 
loop and the Dirac modes, we consider the temporally 
odd-number lattice with the normal non-twisted periodic 
boundary condition for link-variables, in both temporal 
and spatial directions [ 1 ^ [ 2 ^ . 

We introduce the link-variable operator U±^, with the 
matrix element 


(s|17±Js') = 17±^(s)<5s±A..'. (8 ) 


The covariant derivative operator on the lattice is in¬ 
troduced as 

= (9) 

and the Dirac operator 




1 

2 a 




U-b, 


( 10 ) 


with its matrix element 
1 ^ 


^s,s' — ' 7ai U-ij_{s)Ss-fi^s'] : (H) 


fi=l 


where = [/t(s - fl) and 7 I = 7 ^. 

Since the Dirac operator is anti-hermitian, the Dirac 
eigenvalue equation reads 


b\n)=iXn\n), (12) 

where A„ G R. Using the Dirac eigenfunction ipn{s) = 
(s|n), one arrives at the eigenvalue equation 

^ E^=l -k A) - U-f,{s)^pn{s - A)] 

= iXn 1 pn{s)-b^) 


At finite temperature, imposing the temporal anti¬ 
periodicity for D 4 acting on quarks, it is convenient to 
add a minus sign to the matrix element of the temporal 
link-variable operator Lf +4 at the temporal boundary of 
t = Nt{= 0) [il: 

(s, At|174|s, 1) = - 1 / 4 ( 8 , At), 

{s,l|U_4|s,At) = -C/_4(s,l) = -C/l(s,At). (14) 
Then, the Polyakov loop in Eq. 0 is expressed as 


L = — 


^ TTc{Ut}= ^ 


Nt-1 


N,V 


N,V 


^trc{ 1 / 4 ( 5 -knt)}, 
(15) 


n=0 



FIG. 2: The link path structure on a temporally odd-number 
lattice with Nr = 5 and with periodic boundary condition. 
The left configuration is gauge-variant, whereas the middle is 
gauge-invariant. The right configuration represents a closed 
path with At-( link-variables) because of the periodicity in a 
temporal direction, thus this path corresponds to a gauge- 
invariant Polyakov loop. 


where Trc denotes the functional trace, Trc = X)s 
and trc is taken over color index. The minus sign stems 
from the additional minus on 1 / 4 ( 8 , Nt) in Ea. lfHl) . 

Note that the functional trace of a product of link- 
variable operators corresponding to non-closed path is 
exactly zero. Indeed, followed by Eq.®, one obtains 

Trc(f7^il/^, = trc^(s|U^it/^2 

S 

Np —1 Np 

=treJ2^i^iib---U^Npis+ + 

S k^l k^l 

=0, (16) 

with Afc 7^ 0 for airy non-closed path of length Np. 

This is understood from Elitzur’s theorem [s^, that the 
vacuum expectation values of gauge-variant operators are 
zero. 

In the following, we show that the Polyakov loop can be 
explicitly expanded in terms of eigenmodes of the Dirac 
operator. 


A. Relation between the Polyakov loop and Dirac 
modes 

We introduce the following key quantity, 

..s ^ iV —1 

I = Ti-c,^{U4p " ), (17) 

where Trc,-^ = ^^trcti^, and tr.y is taken over spinor 
indexes. From the definition of p in Eq. m, it is clear 

- f.Nr -1 

that U 4 P operator is expressed by a sum of products 
of Nr- (link-variables). 

Note that, on a square lattice, it is not possible to con¬ 
struct any closed loop using the products of odd number 
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of link-variables (see Fig. [2] for illustration). By con¬ 
struction, one considers a square lattice with odd N^., 

thus U 4 P ^ in Eq. (fT71) does not contain any contri¬ 
butions from the products of (link-variables) along 

the closed paths. However, due to the periodic bound¬ 
ary condition in time direction, the only exception is a 
gauge-invariant term proportional to . This term is 
thus related with the Polyakov loop operator. 

Based on the above discussion, and applying Eqs. (II3, 
(fT6)l and ([9|) to Eq. (ITTll . one finds that 


I = 


12V 

(2a)^^-i ’ 


(18) 


thus / is directly proportional to the Polyakov loop. 

On the other hand, since / in Eq. (II3 is defined 
through the functional trace, it can be expressed in the 
basis of Dirac eigenmodes as 

n 

= (19) 


Consequently, from Eqs. m and m one finds that, 
on the temporally odd-number lattice, there is a di¬ 
rect relation between the Polyakov loop and the Dirac 
modes [2i,[2i|, 

i = (20) 


which is a Dirac spectral representation of the Polyakov 
loop. 

The relation (Uni) is a mathematical identity, and is 
exactly satisfied for arbitrary gauge configurations. Con¬ 
sequently, this relation is valid, regardless of whether 
the link-variables are generated in full QCD or quenched 
QCD [ 13 . 

The relation (|2n|l enables to investigate the contribu¬ 
tion of different Dirac-modes to the Polyakov loop. Of 
particular interest is the role of the low-lying eigenvalues 
which are essential to identify chiral symmetry restora¬ 
tion in QCD at finite temperature. 

In principle, it is possible to calculate numerically the 
contribution of the Dirac modes to the Polyakov loop 
through Eq. o by using Monte Carlo simulations. 
However, a very large (4 x TVc x H)^-dimension of a Dirac 
operator implies also a large cost of numerical calcula¬ 
tions. This can be, however, remarkably reduced by use 
of the modified Kogut-Susskind(KS) formalism [^ . 

We rewrite Eq. (1201) into the equivalent form 

i = 25^y:Ay-(n|j/«, (21) 

n 


where the KS Dirac eigenstate jn) is obtained by solving 
the eigenvalue equation, 

r]^Df_,\n) = iX„\n). (22) 


Here, the KS Dirac operator r]^D^ is defined as 




fi=i 


(23) 


with the staggered phase ? 7 ^(s), 

77i(s) = l, rj^{s) = i-iy^+-+^--^ {^i>2). (24) 

Consequently, each Dirac-mode contribution to the 
Polyakov loop is obtained by solving the eigenvalue equa¬ 
tion of the KS Dirac operator whose dimension is now 
{Nc X Vy, instead of (4 x W x Vy, as in the original 
Dirac operator. 

In terms of the KS Dirac eigenfunction Xn{s) = (s|n), 
the KS Dirac matrix element (n\U^\m) is explicitly ex¬ 
pressed as 

(n|I7^|m) = ^(n|s)(s|C/^|s + fl){s + /t|m) 

S 

= '^Xn{syUf,{s)Xm{s + fl). (25) 

S 

From the gauge transformation property of link-variables 
and the KS Dirac eigenfunctioi^ the matrix element 
(n|{7^|m) is gauge-invariant 

Note that the modified KS formalism applied here 
is not an approximation, but is a method for spin- 
diagonalization of the Dirac operator. In this study, we 
do not use specific fermions, such as the KS fermion, but 
apply the modified KS formalism as a prescription to re¬ 
duce the numerical cost. 

The relations of the Polyakov loop and Dirac eigen¬ 
modes in Eqs. (1201) and m are exact. They are valid 
at finite temperature and density, and are independent 
of the particular implementation of fermions on the lat¬ 
tice [ 23 , thus can be used to identify the interplay be¬ 
tween deconfinement and chiral symmetry restoration in 
QCD. 


B. The Polyakov loop fluctuations and Dirac modes 

The expansion of the Polyakov loop in terms of the 
Dirac eigenmodes, formulated in the previous section, 
can be also applied to the fluctuations of the real, imag¬ 
inary and modulus of the Polyakov loop. 

Multiplying Eq. (12(11) by the Z 3 factor one ob¬ 

tains the relation between the Z 3 transformed Polyakov 
loop L and the Dirac modes, 

n 

where k = 0 , ±1 is chosen such that, for each gauge con¬ 
figuration, the L lies in a real sector. 
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Taking the real and the imaginary parts of Eq. (1261) , 
the Dirac spectral representation of the longitudinal and 
transverse Polyakov loops reads 

Lr. = , (27) 

n 

Lt = E (e^""*/'{n|174|n)) , (28) 

n 

respectively, whereas, taking the absolute value of Eq. 
(EOl, the following relation is also obtained; 


Since Eqs. dMl), (EZl), dM]) and dMl) are valid for each 
gauge configuration, the Dirac spectral representation for 
different fluctuations of the Polyakov loop and their ra¬ 
tios are directly obtained by substituting Eqs. d2Zl)-dMD 
to Eqs. dl-dB- As an example we quote an explicit 
expression for the Dirac spectral representation of the 
Ra = XaIxl ratio as 


1^1 


( 2 a) 

12E 


E^"" ^(n|174|n) 


(29) 


J 


Ra = 


^{n\U4\n) 


(E„A^-'Re(e2-W3(n|174|n)))' 


where (x) denotes an average over all gauge configura¬ 
tions. 

The explicit analytic relations of the Dirac spectral de¬ 
composition of the real,imaginary and modulus Polyakov 
loop and their fluctuations are the key results of our stud¬ 
ies. Note here that, like Ea. (l20|) . these relations (pfill - 
(1^ are applicable to both full and quenched QCD, since 
we just use Elitzur’s theorem [33| : only gauge-invariant 
quantities survive. All these relations are derived on 
the temporally odd-number lattice for practical reasons. 
However, a particular choice of the parity for the lattice 
size in time direction does not alter the physics, since in 
the continuum limit, a —>■ 0 and Nt -a oo. any number of 
large Nt should give the same result [^,[^. In fact, sim¬ 
ilar relations are derived also on the even lattice, whereas 
a more compact form can be obtained on the temporally 
odd-number lattice [1^ [2^ . It is however difficult to take 
the continuum extrapolation. For instance, the contin¬ 
uum limit of the Polyakov loop itself is still unsettled 
because of uncertainty of its renormalization. However, 
at least, the ambiguity of the multiplicative renormaliza¬ 
tion of the Polyakov loop can be avoided by considering 
the ratio of the Polyakov-loop susceptibilities 


IV. NUMERICAL RESULTS 


To study numerically the influence of different Dirac 
modes on the Polyakov loop fluctuations and their ra¬ 
tios, we further apply the modified KS formalism. This 
amounts in replacing the diagonal Dirac matrix element 
{n\U 4 \n) in Eqs. (I?7l) - (l^ by the corresponding KS Dirac 


^{n\U4\n) 


A^-^Re (e2-W3(„|U4|„)^ 


(30) 


matrix element (n|? 74 |n) [ 1 ^, as 

(n|?74|n) = 4(n|174|u). (31) 

We analyze the contributions from the low-lying Dirac 
modes to the Polyakov loop fluctuations in the SU(3) lat¬ 
tice QCD through Monte Carlo simulations. In the math¬ 
ematical sense, all the obtained relations (l^- (l^ hold 
for both full and quenched QCD. In this paper, we per¬ 
form SU(3) lattice QCD calculations with the standard 
plaquette action at the quenched level on 10^ x 5 size lat¬ 
tice. Numerical studies are carried out both in confined 
and deconfined phases for different couplings (3 = 
and the corresponding temperatures, T = l/(Wa). We 
use the Linear Algebra PACKage (LAPACK) in diag¬ 
onalizing the KS Dirac operator to obtain the eigenvalues 
Xn and the eigenfunctions Xn(s)- The lattice spacing a 
is determined by the zero-temperature string tension of 
a = 0.89 GeV/fm on a large lattice at each [3 = 6 / 5 ^. 
In fact, we here calculate the static quark-antiquark po¬ 
tential V{r) on 16"'’ lattice at each f3, and fit it by the 
Cornell potential, i.e., the Coulomb plus linear form Q, 
to extract the string tension cr. 

In confined phase, we fix f3 = 5.6 on 10^ x 5 lattice, 
which corresponds to a ~ 0.25 fm and T ~ 160 MeV. We 
also calculate the Creutz ratio x(3, 3) ~ 0.35(2), for an 
estimate of the strin g te nsion or the lattice spacing on 
the precise lattice [ 3 , , in spite of an additional con¬ 

tamination from the Coulomb potential. Nevertheless, 
this value is consistent with the string tension a obtained 
from the potential V{r), and leads to a ~ 0.28 fm. Here, 
the average plaquette value is obtained as ([/^j/)=0.53(2), 
which is consistent with the previous SU(3) lattice stud- 














6 


ies [ 3 ^. In deconfined phase, the simulations are per¬ 
formed at /3 = 6.0 on 10^ x 5 lattice, i.e., for a ~ 0.1 fm 
and T ~ 400 MeV. On this lattice, the average plaquette 
value is (C/^j/)=0.60(2), which is also consistent with the 
previous works [s^. For each value of /3, we use 20 gauge 
configurations, which are taken every 500 sweeps after 
the thermalization of 5000 sweeps. 

Since the Polyakov loop and its different fluctuations 
are expressed as the sums over all Dirac-modes, we divide 
the entire tower of the Dirac eigenvalues into the low- 
and higher-lying modes with the insertion of the infrared 
cutoff A. 

Based on Ea. d^ . we introduce the A-dependent 
Polyakov loops, 


I^Ia 


(2a)^--i 

W 


Hn\U4\n) 

|A„|>A 


for the modulus, and 

iLL)A = Cr ^ 

|A„|>A 

{LT)A = Cr A^-ilm(e2-'=*/3(n|174|n)) . 

|A„|>A 


(32) 


(33) 

(34) 


for the real and the imaginary part, respectively, with 
Cr = {2ai)^--^/W. 

Applying the cutoff dependent Polyakov loops from 
Eqs. (fell . (1551) and (1551) to Eqs. (I5|)-(IS|), we also in¬ 
troduce the A-dependent susceptibilities 


T%x)a 


- {Ya)% 


(35) 


where Y stands for |L|, Ll or Lt, and their ratios 


{Ra)a 


(xa)a 

(xl)a’ 


{Rt)a 


{xt)a 

(xa)a' 


(36) 


To differentiate the importance of the low-lying Dirac 
modes on the properties of the Polyakov loop fluctua¬ 
tions and the chiral condensate, we introduce the cutoff- 
dependent chiral condensate {'ipip)A- In terms of the 
Dirac modes, the ('0 ^/’)a is expressed as [ 2 ^ 


(V’'*/')a 


1 

V 


E 

|A„|>A 


2m 


(37) 


where m is the current quark mass. 

The chiral condensate is strongly affected by the low- 
lying Dirac modes. Taking a typical value for the infrared 
cutoff A ~ 0.4 GeV and the quark mass m ~ 5 MeV, 
leads to a drastic reduction of the chiral condensate 

i?chiral = 0.02, (38) 

in a confined phase at T ~ 0 [U, . 


Figure [5] shows the Dirac eigenvalue distribution 
p(A) = l/V^^(i5(A — A„)) in confined (/3 = 5.6) and 
deconfined (/3 = 6.0) phases. Note here that the near¬ 
zero Dirac-mode density p(A cs 0) is apparently finite in 
confined phase, whereas it is highly suppressed in decon¬ 
fined phase. We show in Fig.|5]the bare chiral condensate 
I {ipfp) I psr S' flavor in the confined phase as a function of 
the quark mass m in the lattice unit. The chiral con¬ 
densate remains finite in the small-m region. Thus, from 
Figs. [5] and m it is clear that the chiral symmetry is defi¬ 
nitely broken in confined phase, whereas it is restored in 
deconfined phase. 






An[a'] 


FIG. 3: The lattice QCD result of the Dirac eigenvalue 
distribution p(A) in the lattice unit: (a) the confinement phase 
at d = 5.6 (i.e., a ~ 0.25 fm) on 10® x5; (b) the deconfinement 
phase at /3 = 6.0 (i.e., a ~ 0.10 fm) on 10® x 5 [ 2 ^. 


In the same spirit, we introduce the ratio, 

Rconf = , (39) 

Ha 

to quantify the sensitivity of the Polyakov loop fluctua¬ 
tions to the particular Dirac modes. When, with some 
A, the ratio stays i?conf — 1, then the low-lying Dirac 
modes below the cutoff A have a negligible contribution 
to the Polyakov loop fluctuations. 

In Fig. [5l we show the Monte Carlo results for i?conf 
in a confined phase at /3 = 5.6, for various values of the 
infrared cutoff A. For the sake of comparison, we also 
show in Fig. [ 5 ] the i?chirai ratio, calculated at the same 
temperature and with the light quark mass, m = 5 MeV. 
The ratios, i?conf and i?chirah indicate the influence of 
removing the low-lying Dirac modes with the infrared 
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m 


FIG. 4: The absolute value of bare chiral condensate 
per a flavor in the confinement phase plotted against the 
quark mass m in the lattice unit. The lattice QCD calcu¬ 
lation is done a.% j3 = 5.6 (i.e., a ~ 0.25 fm) on 10^ x 5. 



FIG. 5: The numerical results for the i?chirai and i?conf ra¬ 
tio from Eqs. (O and (EU. respectively, as a function of an 
infrared cutoff A introduced on Dirac eigenvalues, expressed 
in lattice units. The Monte Carlo calculations have been per¬ 
formed on 10^ X 5 lattice at /3 = 5.6, and for the quark mass 
of m = 5 MeV. 


cutoff A on confinement and chiral symmetry breaking, 
respectively. 

From Fig. [5l it is clear that the i?chirai ratio is 
strongly reduced by removing the low-lying Dirac modes. 
Thus, the low-lying Dirac modes, which are important 
modes for chiral symmetry breaking, are also dominant 
to quantify the chiral condensate. In contrast to i?chirai, 
the i?conf ratio is almost unchanged when removing the 
low-lying Dirac modes even with relatively large cutoff 
A ~ 0.5 GeV. 

Thus, the essential modes for chiral symmetry break¬ 
ing in QCD are not important to quantify the Polyakov 
loop fluctuation ratios, which are sensitive observables to 
confinement properties in QCD. The same result is also 


TABLE I: Numerical results for different Polyakov-loop fluc¬ 
tuations (original), and those without the low-lying Dirac 
modes (IR-removed) below the IR-cutoff A ~ 0.4 GeV. The 
results are obtained in quenched QGD on 10® x 5 lattice at 
P = 5.6 (confined phase) and P = 6.0 (deconfined phase) with 
20 gauge configurations. 


p 


original 

IR-removed 

d=5.6 

T^Xa 

3.475 X 10"^ 

3.470 X 10"'‘ 


T\l 

5.307 X 10"'^ 

5.298 X 10"'‘ 


T\t 

6.005 X lO""* 

5.994 X 10"'‘ 


Ra 

0.6548 

0.6549 


Rt 

1.131 

1.131 

/3=6.0 

T\a 

2.965 X 10"® 

2.965 X 10"® 


T\l 

3.015 X 10"® 

3.015 X 10"® 


T\t 

7.848 X 10"^ 

7.848 X 10"^ 


Ra 

0.9834 

0.9834 


Rt 

0.2603 

0.2603 


found in deconfined phase, as seen in Table HI which sum¬ 
marizes our numerical results on different fluctuations of 
the Polyakov loop and their ratios, obtained on the lat¬ 
tice at /3 = 5.6 and /3 = 6.0, with 20 gauge configurations. 
Note here that the analytical relation (l20l) and the subse¬ 
quent formulae of Eqs. (EZD-dini) hold for each gauge con¬ 
figuration, and the contribution from the low-lying Dirac 
modes to the Polyakov loop L is found to be negligible 
[ 2 ^ . This fact inevitably leads to almost equivalence be¬ 
tween the Polyakov-loop fluctuations and those without 
the low-lying Dirac modes, although the statistical error 
is significant because of the small statistics. 

The differences in the influence of the low-lying Dirac 
modes on the chiral condensate and the Polyakov loop 
fluctuations can be understood semi-analytically. From 
Eqs. dSZli-dMI), it is clear, that the contribution of the 
low-lying Dirac-modes with |A„| ~ 0 is suppressed, rela¬ 
tive to the higher-lying Dirac modes, due to the damping 
factor In fact, a Dirac matrix element {n\U 4 \n) 

does not yield a stronger singularity than I/A^^“^, there¬ 
fore the contribution from the low-lying Dirac modes to 
the Polyakov loop [ 2 ^, as well as, to its fluctuations is 
negligible. Hence, the essential modes for chiral sym¬ 
metry breaking do not contribute to a sensitive probe for 
deconfinement in QCD. Thus, this finding suggests no di¬ 
rect one-to-one correspondence between confinement and 
chiral symmetry breaking in QCD. 


V. SUMMARY AND CONCLUSIONS 

The main objective of these studies was to establish the 
relation between the Polyakov loop and its fluctuations 
with the eigenmodes of a Dirac operator. Based on the 
lattice QCD formalism, we have derived a Dirac spectral 
representation of the real, imaginary and the modulus 














of the Polyakov loop and their fluctuations. Although 
the formulation was done on a temporally odd-number 
lattice, this choice of the parity for the lattice size does 
not alter the physics in the continuum limit with any 
large number of Nr. The analytical decomposition of the 
Polyakov loop and its fluctuations is fully general. It is 
independent from the gauge group, the implementation 
of fermions on the lattice, and is also valid at finite baryon 
density. 

To quantify the influence of Dirac modes over the 
Polyakov loop fluctuations, we have performed Monte 
Carlo simulations in the SU(3) lattice QCD. Our calcula¬ 
tions were carried out with the standard plaquette action 
at the quenched level on (10^ x 5)-size lattice at two dif¬ 
ferent temperatures, corresponding to confined and de- 
confined phases. 

We have shown that the low-lying Dirac modes have 
negligible contribution to the Polyakov loop fluctuations. 
This result is intact both in confined and deconfined 
phases. On the other hand, the low-lying Dirac modes 
are essential, in both phases, to quantify the chiral con¬ 
densate. 

These findings, both in analytical formulas and in nu¬ 
merical calculations, suggest no direct, one-to-one cor¬ 
respondence between confinement and chiral symmetry 
breaking in QCD. However, this does not exclude a coin¬ 
cidence of these two properties in QCD since the abrupt 
change of the ground state from chiral symmetry broken 
to restored phase may drive the onset of deconfinement. 

The above conclusion is based on the numerical sim¬ 
ulations on a rather small-size lattice, being far from a 
continuum limit. Thus, our result on the Polyakov loop 
fluctuations suffers from finite size effects. Such effects 
can certainly modify the values of fluctuations at a given 
temperature, however will not change our conclusion on 
the influence of the low-lying Dirac modes on their prop¬ 
erties. The low-lying Dirac modes have negligible contri¬ 
bution to the Polyakov-loop fluctuations because of the 
damping factor which appears in Eg. (EDI) . Al¬ 

though our numerical calculation was performed at the 
quenched level, the derived analytic formulae are fulfilled 
even in the presence of dynamical quarks. It is one of the 
future prospects to perform full-QCD simulations in the 


present formalism to further justify our conclusion. 

In addition, the derived analytic relations connecting 
the Polyakov loop and Dirac modes are mathematically 
exact for arbitrary odd temporal size Nr- Thus we ex¬ 
pect that our conclusion is rather robust in the contin¬ 
uum limit [ 2 ^, [ 2 ^. Moreover, the ambiguity of the mul¬ 
tiplicative renormalization of the Polyakov loop has been 
avoided in the ratio of the Polyakov loop susceptibilities. 
Yet, it is left as an important but difficult task to extrap¬ 
olate these analytic relations to the continuum. 

In addition to the Polyakov loop fluctuations, there 
are further observables which are linked to deconfinement 
properties in QCD and show abrupt, but smooth change 
across the chiral crossover. One of such observables is 
the kurtosis of the net-quark number fluctuations (STl - fs^ . 
Besides, the QCD monopole, in the maximally Abelian 
gauge, is a relevant degree of freedom in the low-energy 
QCD [ 23 , [2l|, and plays a fundamental role for non- 
perturbative phenomena such as confinement and chi¬ 
ral symmetry breaking. Thus, from the future perspec¬ 
tives, it would be of particular interest to investigate such 
quantities in terms of the Dirac-mode expansion and to 
explore the influence of the low lying eigenmodes on their 
properties near the chiral crossover. 
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